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I. 



INTRODUCTION 



The distortion of a digital transmission, caused by 
noise, may require that techniques be used allowing for the 
detection and eventual correction of errors. In digital 
communications, the most significant performance parameter 
from either a bit or message standpoint is that of probabi- 
lity of error. A means of minimizing the probability of 
error is therefore essential to effective communication. 

Error detection and correction was given great impetus 
by Shannon's paper of 1948, [Ref . 7 ], which extended the 
promise of reliable recovery of digital data perturbed by 
noise (Shannon's Second Theorem). The noisy coding theorem 
provides that messages can be transmitted with arbitrarily 
small error if the source rate is no more than the channel 
capacity. 1 While this sets a goal and the conditions neces- 
sary to attain it, the precise method for obtaining a speci- 
fic scheme was not set forth. However, Shannon did use a 
random coding scheme in proving this theorem. As was shown 
by Rice in 1950, [Ref. 6 ], following Shannon's example, 
choosing codewords randomly leads to the result that as the 
code becomes very long, channel capacity is approached and 



1 Channel capacity is the maximum average information 
input to a channel, which means it is properly matched to 
the channel less the average uncertainty at the receiver 
due to channel noise. 
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the probability of error can be made to decrease exponen- 
tially with code length. But choosing codewords randomly 
is not a practical scheme. 

Until the 1960's the application of channel coding theory 
was slow in its development. The establishment of digital 
circuit technology and the realization of a theoretical 
channel (satellite communication; AWGN) provided a powerful 
stimulus for the utilization of practical error detection and 
correction techniques. The technological improvements have 
continued until today fairly large memories are held by 
small chips, and the components needed to implement powerful 
coding schemes are available in reduced sizes and at reason- 
able prices. 

The objective of this paper is to demonstrate the bene- 
fits derived from a minicomputer (DEC PDP-11/40) simulation 
of channel noise applied to a specific error detection and 
correction scheme. This simulation allows for the determina- 
tion of a probability of bit error for a specifiable noise 
density. Results from the simulation are then used to deter- 
mine statistical trends of a particular error correction 
technique for various internal coding parameters. 

The error detection and correction scheme chosen by the 
author is the Viterbi decoding, [Ref. 3 ], algorithm on a rate 
1/n convolutional code. Both the encoding and decoding 
technique are of current interest in extraterrestrial com- 
munications, with the Viterbi algorithm being principally 
considered for satellite communications. Therefore, in 
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order to provide a complete presentation, the scope of this 
thesis will encompass the basic .principles of convolutional 
encoding and Viterbi decoding, as well as the structure of 
the simulation program, results and conclusions. 

All programming was accomplished on the DEC PDP-11/40 
in machine language, thus presenting the opportunity to 
gain valuable insight into the computer structure while 
completing this thesis. 
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II. BACKGROUND 



Digital communication systems (figure 1) are usually 
designed to minimize the probability of error of the re- 
ceived data (bits) , introduced by noise (gaussian, burst, 
fading, multipath, etc.). Thus, error detection and correc- 
tion may be varied according to the degree and type of noise. 
There are two common forms of error protection: (1) retrans- 

mission and (2) controlled redundancy. The first may be 
applied to systems which are not critically affected by the 
ensuing delay, whereas the second does not require signifi- 
cant delays. This thesis deals with the latter in an 
attempt to gain statistical knowledge of the particular 
coding scheme (controlled redundancy) under study. 

Controlled redundancy techniques are commonly divided 
into two groups: (1) block codes and (2) convolutional 

(tree) codes. Again, owing to the characteristics of the 
channel, one of these two codes may be chosen, along with 
a suitable decoding scheme, to achieve the desired probabi- 
lity of error. A decoding scheme for block codes is batch 
oriented, as is its encoding, using definite algebraic 
operations related to the segmented structure of the code. 

On the other hand, convolutional codes are decoded by a 
statistical procedure due to its continuous (bit by bit) 
nature. 



13 



Source ^ncoSer rTl Modulator ^ ^Transmitter 




X) 

a> 

d 



cd 

* l 

CT co 



w 

rH 

V-I o 
cd xi 

PQ CO 



e 

<u 

4-> 

w 

>> 

CO 

c 

o 

*H 

4-> 

cd 

o 

tH 

c 

d 



o 

o 



cd 

4-> 

00 

•H 

Q 

0) 

C 

o 



M-t 

O 



<D 

XJ 

O 

S 





w 


<u 


<u 


rH 


J-l 


a 

u 


O 

,£> 


5) 


d 


6 


•H 


o 


>> 


& 


CO 


CO 





14 




A. PRINCIPLES OF CONVOLUTIONAL ENCODING 



In 1955, P. Elias first proposed the use of convolutional 
(tree) codes for the discrete time memoryless channel, 

[Ref 2 ] . This technique extended the promise of providing a 
class of codes (linear) whose performance would prove supe- 
rior to that of block codes of the same length, [Ref. 9 ] . 

The development of these codes also gave promise of provid- 
ing for a decoder complexity increasing no more than linearly 
with block length and/or encoder memory. These conjectures 
have been verified, for the most part, by such contributors 
as Fano, Reiffen, Forney and Berlekamp. [Ref. 4] 

Since convolutional codes form a definite discipline 
in coding theory, there are some aspects that will require 
clarification in conjunction with the implementation of the 
code. The next few subsections will accomplish this task. 

1 . Description 

The process of block encoding segments (blocks) of 
an indefinite length input sequence into code blocks pro- 
ceeds in a discrete manner, i.e., the contents of one block 
has no bearing on the encoding of another block. The code 
blocks may be completely changed from the input block or 
they may include the input block with check bits determined 
by a definite algorithm. However, when an input sequence 
is processed through a convolutional encoder the input/ 
output is continuous (sequential bit to bit dependence in 
the encoder) and the output is generated at a specifiable 
number of code bits for each encoder input bit. 
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The term "convolutional" originates from the observation 
that the encoded sequence can be regarded as a convolution 2 
of the input sequence with related generator sequences. An 
example of such a procedure is the binary convolution 
(denoted by *) of two sequences, x and £, where (+) and • are 
binary addition and multiplication, respectively. 

Input sequence: 

x(t Q ) = (x Q x 1 x 2 x 3 x 4 ) = (1 0 1 1 0 ) ; 

Generator sequence: 

5(t 0 ) = (g Q g 1 g 2 g 3 g 4 • • • > = d o o i i . . .) , 

where all bits in x and £ are zero for x_^ and g_^ . 
Then the output (code) bit formed by the convolution of 
x(t Q ) with g(t Q ) is 

Z(t 0 ) = x(t Q )* £(t 0 ) = 0 . 

For time tg + At, the sequences and output are 

x ( tg+At) = (x Q x 1 x 2 x 3 x 4 ...) = (1 0 1 1 0 ...), 

£(t Q +At) = (g Q g 1 g 2 g 3 g 4 ...) = (1 0 0 1 1 ...), 

y ( tg+At ) = (Xg • g Q ) = (1 • 1) = 1 . 

For time tg + 2At, the sequences and output are 



to 



2 The convolution of two functions, f ^ and f 2 , 

' co 

f 1 (t-T)f 2 (T)dT 

/ —00 



is equal 



which means one function is folded in time, then their 
product is integrated as one is shifted pass the other. 
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x(t Q + 2 At) = (x Q x 1 x 2 x 3 x 4 ...) = (1 0 1 1 0 ...), 

£(t Q +2At) — (gpg^g 2^3^ ^ ■*•••■•* = (1 0 0 1 1 •••)» 

y(t 0 + 2 At) = (x Q • g-,^)© • g 0 ) * o . 

And the same procedure is applied to subsequent increments 
of At (one bit shift) . 

This technique extends the concept of block encoding 
to permit memory from bit to bit continuously as against 
memory within the block. Using the operation shown above, 
a conventional convolutional encoder (figure 2) may be 
defined as a linear sequential machine with k-inputs and 
n-outputs, where n > k, usually. This machine is constant, 
linear, causal, and finite state with operations over a 
finite field F, commonly binary. 

The parameter k is the number of input bits which 
enter the encoder in a time increment (At) . These bits 
along with others retained in the encoder memory (finite 
length) , form n-output bits (code bits) . For the remainder 
of this paper k will have a value of 1, but it should 
be remembered that this parameter can take on other values, 
as can be seen in figure 3 . 

The definition specified a finite state machine, where- 
as the example dealt with indefinite length sequences. 

Since an infinite length memory is not practical it becomes 
necessary to limit the basic concept (generator sequence) 
and define some new expressions employed in the actual 
encoder configuration. 
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Figure 2, Conventional (n,k) Convolutional Encoder 
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Figure 3. Representation of the i Input S 
Rate k/n convolutional Encoder (1 



a. Code Rate (k/n) 

The code rate is an expression of the number 
of k-input (1) bits per n-output bits, that is, k/n (1/n) . 
This value may be considered a very close approximation 
to the actual value. The actual value would also include 
a very small number of zero code bits used to terminate the 
code sequence. 

b. Memory Length (M) 

For a practical encoder £ must be limited in 
length. The memory length is a measure of this length and 
is the minimum number of memory cells required to generate 
a code. A representation of an encoder, with a memory 
length of 2, is shown in figure 4. This configuration 
corresponds to the Mealy machine in automata theory. 

c. Constraint Length (K) 

A more conventional expression of the length 
of time which an input bit affects the output sequence is 
that of constraint length. This term is simply the memory 
length plus one, and is depicted as the number of memory 
elements in a Moore machine configuration of an encoder 
(figure 5) . As may be expected this value is one factor 
in determining the complexity of the error detection and 
correction scheme. 

d. Free Distance (d^) 

The term distance, as applied to coding, refers 
to the number of differing bits between two code sequences 
of the same length, as seen below: 
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Binary Convolutional Encoder. 



code sequence 1: (1001011001) 

code sequence 2: (1011001001) 

Differences between the two sequences exist in the third 
and sixth bits only, therefore the distance between the two 
sequences is 2. When using block codes to construct a code 
sequence the minimum distance between all pairs of codewords 
(code blocks) is a definite indication of the error detect- 
ing and correcting capability of the code. However, if the 
code does not employ a block configuration, the use of this 
value may not be justified to indicate the codes error 
detecting and correcting capabilities. 

This is the case when discussing convolutional codes. 
Free distance is the term for the value expressing the 
theoretical error correcting capability of a convolutional 
code. This value is defined as that minimum number of set 
bits occurring in a code sequence which resulted from the 
input of a nonzero sequence. Table I shows how the free 
distance (d^) can be determined for a rate 1/2 convolu- 
tional code with K=3. As can be seen, the larger the 
value of K then the number of possible input sequence to 
be considered also increases at a great rate. But with the 
aid of computers the tedious grinding process is accomplished 
rather quickly. 

The convolutional encoder is constructed to take full 
advantage of the free distance applied to a specific code 
rate and constraint length. This is accomplished through 
the proper selection of generator sequences. However, it 
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should be noted that with the computers currently available 
the drudgery can be passed on to the machine. 

2. Properties 

The properties of the encoder, as related to the 
generator sequence, are discussed in this section. For 
convenience, the set of generator sequences (cj^ , . . . , c[ n ) 
used to generate the convolutional code is denoted by 
[G] , a matrix representation of the encoder discussed in 
the next section. Figure 6 is useful in relating the 
following discussion to a physical communication system. 

a. Property 1 

Foremost, the useful encoder should generate 
a code which will yield the fewest number of errors in the 
codeword estimator. The estimator is a demodulation scheme 
at the receiver determining as accurately as possible the 
received sequence before actual decoding takes place. Along 
with this error minimization, the complexity of the estima- 
tor requires minimization to develop a useful error correc- 
tion system. Error minimization and complexity minimization 
together may require some compromise in the practical 
system depending on the desired reception quality and sys- 
tem environment. 

b. Property 2 

The encoder [G] must have an inverse relation 
for decoding purposes. This is realized by the relation 

y[G] _1 = x [G] [G] -1 = D P x 
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is physically 



for all x. The matrix (decoder) , [G] ^ 

realizable and is a pseudo inverse, in that, the expression 
[G] [G] ■*■ , yields a delay D^, where p is the number 
of time intervals before decoding occurs in the expression. 
If p is zero the decoding occurs on receipt of the 
sequence and the above relation becomes 

^[G]" 1 = x [G] [G]” 1 = x, 

where [G] ^ = [G] ^ when p = 0. 

A mention of the possibility of catastrophic error 
propagation, [Ref ] , is required at this time to stress 
the point that [G] ^ must be feedback-free. This may be 
interpreted as meaning that the n generator sequences, 
in polynomial form should not have a common factor. If a 
common factor exists then the input of a sequence with a 
finite number of set bits to the encoder, may be decoded, 
after noise is added, as a sequence with an. infinite 
number of set bits. This is catastrophic error propagation 
and is avoided by making [G] , thus [G] ^ feedback-free, 
c . Property 3 

The matrix [G] must also meet the requirements 
mentioned in the definition of a convolutional encoder. 

(1) Constant (time invariant) . The time 
invariance of the encoder is represented by 

[G] (D P x) = D P y , 
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which means if all the inputs are shifted in time, then 
all the outputs are shifted accordingly. 

(2) Linear. The output sequence resulting 
from the superposition of two input sequences, x, must 
be equal to the superposition of the two output (code) 
sequences, y_, that would result from the inputs entered 
separately. This is also necessary for the multiplication 
of an input by a scalar as shown below. If G:x->^, then 
G(x 1 +x 2 ) = GCx^) + G(x_ 2 ) = Y-i + Y.2' and G ( ax^ ) = aGfx^) = 
a Xl w here a is an element in the field F (usually GF(2)). 

(3) Causal . The existence of a nonzero output, 
y, prior to the input of a nonzero sequence, x, is for- 
bidden. This is accomplished by the encoder output being 
zero when its memory elements are all zero, therefore zero 
inputs. 

(4) Finite State . The states of the encoder 

are finite in number due to the value of the constraint 

K — 1 

length. Thus, a binary encoder has 2 possible states. 

Each state being those input elements involved in the 
generation of the next output (code) bit, when combined 
with the next input bit entering the encoder. 

d. Property 4 

The chosen encoder [G] should have the minimum 
number of memory elements to generate the code. 

e. Property 5 

The chosen encoder [G] should generate a code 
yielding the fewest number of decoding errors per error 



28 



event (channel noise) . Thereby, the fastest correct 
decoding decision is made at the receiver. 

The preceding paragraphs provide a basis by which the 
acceptable class of encoders is greatly reduced in number. 
However, the selection of [G] in its optimum form for a 
given channel is a long drawn out process, possibly best 
suited for computer analysis. Theoretical selections have 
been made by many people, [Ref. ] , who based their deci- 
sions on the free distance alone. 

3 . Representation 

The encoder configurations used to this point are 
accurate physical descriptions, but they lack a convenient 
form needed for analysis. There are five encoder repre- 
sentations which show varying degrees of the code's struc- 
ture. Each is discussed in the following sections. 

a. Polynomial 

Using the D-transform, the input and output 
sequences appear as polynomials of degree much greater then 
the constraint length (K) . 

-1 2 

Input, X(D) = ... + x_^D + Xq + x^D + + • • • • 

output, (d) = ... + y^^D 1 + y i/0 + y i,l D 

2 

+ y . „D + . . . , 

where 1 £ i .£ n. The generator sequence, for the Mealy 
machine configuration in figure 7, is a polynomial of 
degree equal to K-l. The mapping of X -*■ is accom- 

plished by polynomial multiplication. 
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Figure 7. Polynomial Representation of a Generalized 
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where L is the bit length of the input sequence. This 
representation has little advantage in a code structure 
analysis, but it is very definitive as to the procedure for 
code generation. 

b. Matrix 

A matrix, [G] representation of a rate 1/n 
encoder starting at state zero and time zero is achieved 
by the manipulation of each generator sequence. When 
y = [G] x , [G] may be shown as a matrix whose row vectors 
are the generator sequences (figure 8a) . Using this nota- 
tion restricts x to be a length K. However, the require- 
ment to change x for each code computation is eliminated 
by the matrix in figure 8b. This notation incorporates 
the shifting of x past each £^(1 ^ i ^ n) into the 
matrix [G] . For example, (figure 8 follows the example). 
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n = 2 
K = 3 

£]= (1 0 1 ) 
£ 2 = (1 1 1) 



[G] 



11 01 11 00 00 . . . 

00 11 01 11 00 . . . 

00 00 11 01 11 . . . 

00 00 00 11 01 

. 00 11 
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This form allows for a fast generation of the code and 
is definitely a model which could be implemented on a com- 
puter with little difficulty. Again, there is a lack of 
insight into the code structure except for the possible 
application of matrix algebra theory, 
c. Tree 

The most common representation is a tree dia- 
gram (figure 9) which incorporates the branch and nodal 
properties of the code generation. The base branch corres- 
ponds to the initial state of the encoder prior to a non- 
zero entry. The first node (a) is state zero for the en- 
coder. When a 1 begins the sequence the first lower 
branch is chosen and the code is found on this branch 
leading to node (b) or state 1. The subsequent input 
bits dictate whether the tree is followed up or down. 
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Figure 8. Generalized Matrix Representation of a Rate 
1/2 Convolutional Encoder (K=3) . 
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corresponding to inputs of 0 and 1 respectively. If 
this is continued to the branch, all of the possible 

input sequences of K bits have been represented, thus 
each possible encoder state. Beyond the branch it is 

readily recognized that the code symbols on the branches 
leaving from the two nodes labeled ( a) are identical iden- 
tical to those at the base of the tree. Even though the 
structure has become recurrent, for long sequences this 
diagram may become too large to analyze easily, 
d. Trellis 

Since the input sequence (1 0 0 ...) and 
(0 0 0 ...) generate the same n-bit codewords after the 
third branch (figure 9) , then both nodes labeled a can be 
joined together. The recursive character of the tree 
diagram lends itself to be redrawn with remerging branches, 
thus forming a trellis (criss-cross) diagram (figure 10) . 
The encoder state is the basis for the trellis diagram. 
These states represent those input bits in the memory cells 
which will generate the code bits when the next input bit 

K_ ]_ 

enters the encoder. Therefore, there are 2 states on 

either side of the diagram representing the present and 
next states. The branches connecting the states denote the 
code bits being generated when a 1 (dashed line) or a 
0 (solid line) enters the encoder. Again, after the third 
input bit the complete trellis diagram is specified and is 
repeated for all succeeding branches. 
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Figure 9. Tree Representation of Convolutional 
Encoder in Figure 5. 
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Figure 



For short constraint lengths (K < 10) , this repre- 
sentation is highly desirable due to its compactness. The 
trellis diagram is used in the Viterbi decoding algorithm 
and easily lends itself to be listed in a computer program, 
e. State Diagram 

The last encoder representation is similar to 

the trellis, but is spread out more to allow more parameters 

to be placed on the branches between states (nodes) . The 

state diagram (figure 11) has the input 0 and 1 denoted 

in the same manner as for the trellis, but some additional 
* 

notation is placed on' the branches. 

(1) D^. The value of q denotes the weight 

of the code bits for that branch. 

(2) L^. The value of r denotes the branch 

length of a path from state x to state y. 

(3) N S . The value of s denotes the number 

of input one branches encountered in a path from state x 
to state y. 

All of the above notation is used in figure 11 and the 
results of a state path for the input sequence (10100) 
are shown . 

Since the state diagram is a directed graph, a transfer 
function can be determined using the theory of directed 
graphs, [Ref. 4, pp 2A-1-2A-11] . The transfer function con- 
sists of various powers of the three measurements listed 
above. These values are used to determine the properties 
(D, L, N) for all paths in the convolutional code, when the 
transfer function is in expanded form. 
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Figure 11. State Diagram Representation of 

Convolutional Encoder in Figure 5. 
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The trellis and state diagrams are the most complete 
representations and offer the analyst a chance to study 
the path structure of the convolutional code. Each diagram 
yield a compact representation, easily arranged for small 
values of K. However, the looping branches necessary for 
the construction of a state diagram become very confusing 
and difficult to arrange for larger constraint lengths. 
Therefore, the trellis diagram is not only easy to construct 
and understand, but if the values of D, L, N, are needed 
for analysis they can incorporated into the diagram in an 
orderly fashion. 

The discussion of convolutional encoding principles 
was basic and is by no means a complete detailed course in 
convolutional codes. However, the purpose of this section 
is to introduce the reader to concepts which will aid him 
in his understanding of the thesis computer program. 

B. PRINCIPLES OF VITERBI DECODING 

Decoding is the inverse operation to encoding and is 
intended to recover the source bits with all, or almost 
all, of the channel errors removed. The decoder, which 
may be implemented as hardware or software, utilizes the 
encoded bits to detect and/or correct errors. Error detec- 
tion is similar in complexity to the encoding operation. 
Error correction, however, is inevitably a more complicated 
process than encoding, since the goal is to reduce the 
probability of a decoded bit error. 
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A scheme for decoding convolutional codes was proposed 
by Viterbi in 1967, [Ref. 8]. Tt was shown by Viterbi that 
an optimum decoding procedure existed for a statistically 
independent (bit to bit) input sequence transmitted over a 
channel whose errors occur independently from channel bit 
to channel bit. Subsequently, Forney found that the Viterbi 
decoding algorithm is synonomous with maximum likelihood 
sequence decoding, [Ref. 3]. 

1 . Description 

Simply stated, the Viterbi algorithm is a solution 
to the problem of finding the most likely encoded sequence 
through a state (finite) diagram representing the encoder. 
For a decoder to minimize the overall error probability 
of a decoded bit by brute force, maximum likelihood decoding 
would mean calculating the likelihood of the received 
sequence on all paths of the encoder state diagram. How- 
ever, there are two factors which reduce the complexity of 
this problem. The first is the fixed periodic structure 
of the encoder trellis diagram, and the second is the code 
characteristic of remerging paths after the same K input 
bits are applied to two different paths (figure 12) . For 
these reasons, the trellis diagram is an ideal tabular 
representation of the flow of the code at any instant. 

The decoding process, as mentioned above, was found 
optimal for a statistically independent input sequence in 
discrete time. Along with this stipulation, the channel 
noise is also memoryless, as is the case for a binary 



41 



Paths 

Separate 




w 


/-N 


•H 


CO 


r— I 


!! 


i — 1 




<U 




U 




H 


U 


C 




•H 


XI 




o 




a 


4-» 


c 


cd 


w 


Ph 




60 


cd 


c 


c 




o 


60 


•H 




4-» 


0) 


31 


e 


i — 1 


0) 


O 


P4 


> 




c 


4-1 


o 


O 


CJ 


0) 


4-1 


1 — 1 


O 


a 




B 


fi 


cd 


cd 


K 


M 


w 


60 




cd 


s 


•H 




Q 


CM 




1 — 1 




0) 




U 




3 




60 




*H 









42 



symmetric channel (BSC) with hard decision demodulation 
(2 1 levels) or an additive white gaussian noise (AW3N) 
channel with soft decision demodulation ( 2 ® levels, Q > 1) . 
The BSC errors transform a 0 to a 1 and a 1 to a 0 
and occur independently from bit to bit with probability 
p. In the AWGN channel, the probability of a given quan- 
tized value (0 to 2*^-1) of a received bit is determined 
from the gaussian probability density function. These 
values of the received bits are used to determine the most 
likely received sequence. 

A scoring procedure is employed indicating the trellis 
path which shows the least difference from the received 
sequence. This is accomplished by using the concept of 
distance for the BSC and a numerical difference for the 
AWGN channel. The latter may be called the innerproduct 
of the received bit and the calculated (trellis) trans- 
mitted bit. This is implemented digitally by taking the 
exclusive OR of the two Q bit representation of the bits 
above. Both (BSC and AWGN) techniques yield a value that 
can be used to determine the most likely (lowest score) 
path . 

2 . Implementation 

In order to describe effectively the implementation 
of the Viterbi decoding algorithm an example follows with 
step by step explanations supplemented by appropriate dia- 
grams. This example will be for a rate 1/2 convolutional 
code with K = 3. The generator sequences are 5 Q (1 0 1) 



43 



and 7g (1 1 1) . The following sequence is the encoder in- 
put: (1 01100011Q100010 ...). 

The corresponding code sequence is 
(11 01 00 10 10 11 00 11 10 10 00 01 11 00 11 01 11 00 ...). 
An error sequence denoting BSC noise perturbs this sequence 
prior to being received at the decoder. The error sequence 
is (01 00 00 00 01 00 00 00 00 10 00 00 00 00 00 10 00 00 



...). Then, the received sequence is (10 01 00 10 11 11 00 
11 10 00 00 01 11 00 11 11 11 00 ...). 



The following subsections are the steps of the Viterbi 
algorithm decoding the perturbed (received) sequence, 
a. Step 1 

The following figure and tables are set up to 
indicate the trellis structure, the path (survivor) sequences 
and the path scores. 

(1) Trellis . (Figure 10). 

(2) SSEQ (I) . The present state survivor 



sequence. 

(3) SCORE ( I ) . The present state survivor 
sequence score. 

(4) SSEQ ( J) . The next state survivor sequence. 

(5) SCORE (J) . The next state survivor sequence 

score. 

b. Step 2 

An initialization of tables (2) and (3) is 
shown in figure 13. Tables (4) and (5) are not considered 
until the first group of 2 (n) received bits enter the 
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decoder. The all zero sequences placed in SSEQ(I) depicts 
the assumption that the encoder had no set bits entered 
prior to the beginning of the message under consideration. 

A decoder parameter is the length of the survivor sequences 
used to hold path estimates until a bit is decoded with the 
desired probability of error. This length is called the 
decoder constraint length (DCL) , whose value is 6(2 x K) for 
this example. Also initialized is SCORE (I) so that a score 
of zero is in state 0 to denote that the encoder is 
assumed to have started from this state. The other states 
are assigned scores which demonstrate the unlikelihood of 
the encoder starting from these. Normally, a value of 
2 x K is sufficient, which is 6 in this example, 
c. Step 3 (Figure 14) 

This step is the first in the actual decoding 
process. However, it will become obvious that it is recur- 
sive throughout the remainder of the algorithm. Now, the 
first 2 received bits (10) enter the decoder. These bits 
are compared with each pair of code bits on the branches 
of the trellis and a distance (A) is determined for each 
branch. This A is then added to the present score of the 
corresponding present state. Now there are two scores 
coming in on the branches to each next state. At this time, 
a decision is made as to which branch, coming into each next 
state, has the lowest score. If the scores are identical 
then an arbitrary choice can be made, decode estimate is 
zero for this example or an alternative may be exercised to 
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further evaluate the A added to the score in SCORE (I). 

When all the next states have a score, they are entered in 
SCORE (J). Then the next state's survivor sequences are 
determined for entry into SSEQ(J) . This is accomplished by 
determining what bit, 0 or lwas dropped from the present 
state at the end of the minimum score branch for each next 
state. Then, this bit is shifted into the beginning of the 
appropriate present state sequence (SSEQ(I)) and the entire 
sequence is then placed in the SSEQ(J) table, corresponding 
to its next state. 

The oldest bit just shifted out of the sequence in 
SSEQ(I) at state 0 is the most likely decoded bit since 
it is associated with the minimum score (0) sequence in the 
present state. The bit decoding, after a length of DCL, 
minimizes the effect of randomly spread errors in the 
received sequence. Again, if the scores of two or more 
sequences are equal, then a decision is made according to 
a more detailed comparison of the scores in SCORE (I). 
d. Subsequent Steps 

Since the most recent paths and scores are now 
in SSEQ(J) and SCORE (J) , the roles of present and next are 
swapped, (I)<— > (J) , and step 3 is carried out again. Thus, 
after each decoding bit decision, the (I) and (J) are 
swapped, and step 3 is repeated time and again. Figure 15 
shows the decoding steps for four more shifts of the decoder 
input. The tables, after the eighteenth shift, are shown 
in figure 16 . 
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Figure 15. The Next State Tables for the Second Through Fifth 
Shift into the Decoder. 
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It is noticed, at this time, that the first 8 
(DCL + K-l) bits are zero. This is explained by the pre- 
sense of 6 (DCL) zeroes in the initial SSEQ(I) and the 
2 (K-l) zeroes in the encoder when the first nonzero code 
bits were generated. 

Another note should be made of the ease with which a 
soft decision demodulation scheme could be incorporated 
into the scoring process. This dimension would produce a 
more accurate representation of the likelihood of a sequence 
and possibly eliminate the need for arbitrary decisions, 
which may also be time-consuming. 

At a first glance this procedure may seem strange or 
awkward, but 'there can be no doubt that the Viterbi algo- 
rithm is easily computer (hardware or software) implement- 
able. The recursive nature of the primary decoding step 
is the prime factor in controlling the size of that imple- 
mentation. 
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III. COMPUTER PROGRAM 



The term simulation in the title of this thesis should 
not be construed to mean the encoder and decoder operations 
(Programs) are simulated. By no means is this the case, 
for the rate 1/n convolutional encoder and the Viterbi 
decoder are software implementations that could be used 
in an actual system. Of course, some form of synchronization 
is needed for practical operation, but in this program that 
is assumed to have been accomplished by elements (hardware 
or software) preceding the decoder in the receiver. 

A. CHANNEL NOISE 

Actually, the simulation occurs when a channel is des- 
cribed in the program by a noise generating section. This 
section includes a quantization segment, random number seg- 
ment, noise generation segment, and perturbation (summing) 
segment. The following paragraphs are discussions of these 
four divisions of the channel noise program. 

1 . Quantization 

This part of the program uses the encoded sequence 
as an input, then passes this sequence on to the noise sec- 
tions after Q zeroes or ones have been substituted for 
each 0 or 1. When Q=1 is a program input this segment 
is bypassed for apparent reasons. 

2 . Random Numbers 

The random numbers generated in this program were 
obtained by using the Lehmer congruential method. 
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X n+1 = aX n + b C” 0 - 3 T o» ' 

where a = 257, b = 1, and Tq = 2^. The term Xq, starting 

number, is varied to provide representative sequences of the 

distribution shown in figure 18. These numbers are used to 

determine the time between set bits in noise sequence. The 

T 0 

period of the random number sequence is 2 , therefore, 

T 0 

a large sampling of the sequence approaching 2 would 
result in a binomial distribution. However, the test runs 
used to determine results in this thesis employed message 
inputs of 10,000 bits in length. Thus, the distribution 
(dots) shown in figure 17 is a better representation of the 
time between set bits in the noise sequence. 

3 . Noise Generation 

After a random sequence has been generated in the 
previous segment, 2. , a noise sequence is formed using 
the numbers. Multiplication of the number of set bits in 
16 bits (T) , the length of one random number, with powers 
of 2(2 2,2 ^,2^,2^) determine how many zeroes will occur 
before the next set bit in the noise sequence. For example, 
if the random number being considered by the program is 
0110011100100100, then its weight is 7g or lll^. If the 
noise parameters (core locations 10230 and 10236) specify 
multiply by 2 ^ then the number 7g is shifted one place 
right and becomes 3g. Likewise, if multiply by 2 1 is 
entered, then 7 g becomes 16 g . These new numbers are then 
used to determine the space between set bits in the noise 
sequence. Each 16-bit number from the congruent ial random 
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Figure 17. Comparison of Congruential Random Number Generator 
(10"* bits) with Poisson Distribution. 



number generator is used to generate a space. The number 
of 16-bit numbers used in one run is determined by the num- 
ber of code bits which are present when the noise is added 
to the code sequence. 

The variation in the set bit density for the noise 
sequence is analogous to varying the density of channel 
errors. This fact is used in the evaluation of the codes 
in section IV. Figures 18 through 23 show the probability 
of k errors in 16 channel bits (Q = 1,2,3) versus the 
number of errors, k, in 16 channel bits. The plots are 
for averages (XT) of about 10 and 2 errors per 16 channel 
bits. Also shown is the corresponding plot of the Poisson 
distribution as a continuous curve. Any discrepancy between 
the theoretical and measured (dots) is attributed to the 16 
bit length constriction primarily. 

4 . Perturbation 

The purpose of this last segment of the noise 
simulation is to simply exclusive OR (additive noise) the 
noise sequence to the quantized encoded sequence. This 
operation forms a perturbed sequence, which is the received 
sequence for the decoder. If this was to be changed to 
multiplication or some other operation another noise form 
could be simulated for the coding scheme. 

All of the details, as to programming specifics, are 
discussed in Appendix B, along with a program listing. 
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Figure 18. Simulated Noise ( Q = 1, AT = 9.8) 

vs Poisson Distribution (AT = 9.8). 
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Figure 19. Simulated Noise (Q = 1, XT = 2.0) 

vs Poisson Distribution (XT =2.0) 
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Figure 20. Simulated Noise (Q = 2, AT = 9.8) 
vs Poisson Distribution (AT = 9.8 
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Figure 21. Simulated Noise (Q =2, XT = 2.0) vs 
Poisson Distribution (XT = 2.0) 
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Figure 22. Simulated Noise (Q = 3, IT = 9.8) 
vs Poisson Distribution (XT = 9.8) 
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Figure 23. Simulated Noise (Q = 3, XT = 2.1) vs 
Poisson Distribution (XT = 2.1). 
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B . PROGRAM FLOW 



The four major sections of the program are the encoder, 
decoder, noise generator, and the message generator 
(ASCII 7-bit code) . An appendix is devoted to a detailed 
listing and brief discussion of each of the above sections. 

The parameters that are needed to initiate the simula- 
tion are: 

1. N; the inverse rate (1/N (memory location 10200 
in Appendix F) . 

2. K; the encoder constraint length (memory location 
10206 in Appendix F0. 

3. Q; the quantization levels (memory location 10214 in 
Appendix F) . 

4. DCL; the decoder constraint length (memory location 
10222 in Appendix F) . 

5. Noise parameters; operation (-j-/x) and operand 
(1 , 2 , 3 , . . . ) (memory locations 10236 and 10230 
respectively in Appendix F) . 

6. Generator sequences; representations entered so 
that first cell in encoder is bit 0 (memory 
locations 13760-13776, where G ± (13760, G 2 (13762, ... , 
Gg (13776) ) . 

Another programming concept, position independent code 
(PIC) , is used to allow the user to move the entire program 
or any part (with few changes) to another part of core. 

This approach introduced the need for a program segment that 
initializes all addresses and counters (constants) used in 
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the program. The listing and a brief description of this 
segment are included in Appendix D. 

The entire program, in block form, is depicted in 
figure 24. If the reader is interested in a detailed 
program description he is referred to Appendices B-E. 
Appendix A contains a map of that part of core used for the 
implementation of the position independent code. 
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(10000-10170) 



Message Source 
ASC II 7 Bit Code 



(10174-10706)i 



Address and Counter 
Initialization 



(10710-11176)^ 




Rate l/n Convolutional 
Encoder 


(11200-11706x1 


l 


Noise Simulator 


(11710-1317 6) J 


r 


Viterbi 


Decoder 


(13700-13400) V L 


Message 
ASC II 7 


: Sink 
Bit Code 



Note: The numbers in parentheses are 

actual addresses occupied by 
the instructions corresponding 
to a particular flow chart 
block. 



Figure 24. Block Diagram of Program Flow. 
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IV. RESULTS AND CONCLUSIONS 



The term best code is used in many papers, but for the 
most part the use of the term best code is not defined in 
sufficient detail. This is one descrepancy that will be 
avoided by this thesis. Included in this section is a 
summary of the computer results obtained from a best code 
(defined below) determination procedure and encoder/decoder 
parameter variations. 



A. BEST CODE DETERMINATION 

The determination of a best code for a specific channel 
a simulated channel in this paper, follows a procedure based 
on Shannon's fundamental theorem for a discrete channel with 
noise. This procedure is best described using Shannon's 
representation of the attainable region in a graph (figure 
25) of H(x) (information rate) versus H^(x) (probability 
of decoded messager error). If H(x) < C (channel capacity), 
then Shannon shows that H (x) can be made arbitrarily small 

y 

with a proper encoding procedure. When H(x) > C, then the 
excess information being pushed onto the channel can only 
increase the uncertainty (H^.(x)) of the decoded message. 

The minimum value of H^(x) is very close to H-C in this 
latter case. 

The computer program in this thesis is implemented with 
a supplemental segment which stepped through all well-defined 
convolutional code generators of a specific code rate and 
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constraint length. Well defined refers to the constraint 
length definition of the code, meaning that the coefficient 
of the zero power in the generator polynomial must be set for 
at least one generator sequence, and then another sequence 
must have the coefficient set for the K-l power. An 
example of this is shown here: 



n 



n 



2, K 



3, K 



4; 2 1 (D) 

2 2 (D) 
3; 

22(D) 

2.3 (°) 



2 3 
D + D + D 

3 2 
D + D =*> 1 + D 

2 3 2 

D + D =$> D + D 

D + D 2 + D 3 =£> 1 + D + D 2 



The rate is defined by the number of nonzero generator 
sequences . 

The first two paragraphs describe a procedure and the 
selection of encoders to be processed by the program. The 
program is applied to find the generator sequences, for a 
given code rate/constraint length, which come the closest to 
the lower error (H^.(x)) boundary of figure 25. Code rates 
of 1/2 and 1/3 are tested and if the input message is 
assumed to have 1 bit of information in each bit entering 
the encoder, then these rates correspond to those plotted 
on the horizontal scale of figure 25. Tables II through IV 
present a summary of the test run results for a best code 
determination . 

Each test run was made for a 40,000 bit encoder input 
applied over various samples of the simulated noise and at 
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Measure of average error 



A 




Information Input Rate 



C = Channel Capacity 
H m = Maximized H (x) 
ax at Transmitter 

Figure 25. The measure of average error for a given 
information input rate to a channel 
(figure is similar to Shannon's 
Representation) . 
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Encoder Constraint Length (K)=3 
Decoder Constraint Length (DCL)=16 
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four values of XT(1,2,4,10 errors per 16 bits). The com- 
puter output is a listing of the decoded message errors for 
each well defined encoder. This output would require many 
additional pages to list in this thesis. Therefore, the 
summary in Tables II through IV is presented instead. The 
number of decoded errors is then divided by 40,000 to give 
the measure of average error for each test run (AT) . 

The various values of XT are obtained by the noise 
simulation and have the characteristic distributions of those 
curves presented in figures 18 and 19. 

A discussion of these results is presented in part C of 
this section. 



B. DECODER PARAMETER VARIATIONS 

In order to provide an insight as to the effects of 
quantization (Q) and decoder constraint length (DCL) , tables 
V and VI are shown. The effect of increasing Q or DCL is to 
decrease the frequency of errors thus improving the error 
correcting scheme. The codes used in these tables are those 
chosen as the best codes from tables II through IV. 



C. DISCUSSION AND CONCLUSIONS 

The determination of a best coding scheme for a specific 
channel application can be made from tables II-IV. The 
generator sequences which have the minimum frequency of 
errors (in place of (x) ) for the various values of XT 
are outlined boldly. Using this comparison technique plus 
a comparison of the general statistics of the encoder at 
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other values of XT, the following generator sequences are 
chosen to be the encoders which meet the criteria set forth 
previously. 



<N 

II 

S 


K = 3; 


% = <H0) . £2 = 


(001) 


N = 2, 


K = 4; 


Sl = (mo) , a 2 = 


(0101) 


N = 3, 


K = 3; 


£l = (in) , a 2 = 


(001) , 






3-3 m 


(001) 



These encoders were also used for the results obtained by 
varying Q and DCL in tables V and VI. 

With actual concrete measurements of the effectiveness 
(minimum error) of a code and a means (computer program) of 
applying these codes to actual channel recordings there is 
no reason to have to settle for an inadequate error correcting 
system. The results presented in this section indicate clear 
preferences in choosing certain codes to accomplish desired 
communications in a given noise environment. The user speci- 
fics the task and with the computer a code can be chosen. 

The true value of such a program should be evaluated when 
implemented with actual channel noise. No one has been able 
to find a relation or algorithm for encoding that would enable 
a communication system to reach the maximum average informa- 
tion rate (c) with an arbitrarily small frequency of error. 
However, computer analysis offers the channel user an oppor- 
tunity to find a coding scheme which meets the requirements 
for that channel. 
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Decoder Constraint Length = 16 Bits 
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APF’ENDI X fi 



POSITION INDEPENDENT CODE CORE MAP 

THE CORE LOCATIONS OF ADDRESSES, COUNTERS, AND 
MISCELLANEOUS STORAGE (WORK) AREAS ARE LISTED BELOW. 
SOME OF THESE VALUES, PRECEDED BV *, ARE ENTERED BV 
OR COMPUTED BV THE PROGRAM INITIALIZATION SUBPROGRAM, 
WHICH FOLLOWS THE INPUT SUBPROGRAM. 



* 13600; 


INPUT BLOCK ADDRESS 


* 13602; 


ENCODED/DECODED BLOCK ADDRESS 


* 13664; 


QUANTIZED CODE SEQUENCE BLOCK ADDRESS 


* 13606; 


NOISE SEQUENCE BLOCK ADDRESS 


* 13610; 


RANDOM NUMBER/ PERTURBED SEQUENCE BLOCK ADDRES 


13612; 


TRELLIS TABLE BLOCK ADDRESS 


* 13614 ; 


SSEQ(I) TABLE BLOCK ADDRESS 


* 13616; 


SCORE* I > TABLE BLOCK ADDRESS 


* 13620; 


SSEQ(J) TABLE BLOCK ADDRESS 


* 13622; 


S C 0 R E ( J ) TABLE BLOCK ADDRESS 


13624; 


"NOT USED" 


* 136 2 6 ; 


GENERATOR SEQUENCE BLOCK ADDRESS 


* 13630; 


STACK ADDRESS 


13632; 


HOLD ( STORAGE) 


13634; 


DELTA (USED IN DECODtR) 


1 3 6 3 6 ; 


MASK (USED TO MAKE DECODED BIT DECISION) 


13640; 


WORK (STORAGE) 


13642; 


MESSAGE WORD (16 BITS) COUNT WORK LOCATION 



13644; MG! BIT COUNT WORK LOCATION 
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12646; 


DECODED BIT COUNTER 


12650; 


NUMBER OF SCORES COUNTER 


1 2 6 5 2 ; 


MINIMUM SCORES ADDRESS 


12654; 


MINIMUM SCORE VALUE 


12656; 


N-BIT LOOP COUNTER 


12660; 


ZERO BIT SCORE VALUE 


12662; 


ONE BIT SCORE VALUE 


12664; 


SSEQ WORK LOCATION 


12666; 


TOTAL OF NOISE VALUES IN NOISE SEQUENCE 


12670; 


TOTAL NUMBER OF ERRORS IN NOISE SEQUENCE 


12672; 


TOTAL NUMBER OF DECODED MESSAGE BIT ERRORS 


12674; 


"NOT USED" 


12676; 


"NOT USED" 


* 12780; 


NUMBER OF CODE BITS PER INPUT BIT, N 


* 12702; 


ENCODER CONSTRAINT LENGTH, K 


* 12704; 


QUANTIZATION VALUE, Q 


* 12706; 


DECODER CONSTRAINT LENGTH, DCL 


12710; 


MESSAGE WORD <16 BITS) COUNT, W 


* 12712; 


POWER OF TWO USED IN GENERATING NOISE SEQUENCE 


* 12714; 


OPERATION <DI VI DE/MULTI F'LV) USED FOR NOISE 


* 12716; 


Q X N 


* 12720; 


Q X N X W 


* 12722; 


4 X K < NONZERO STATES INITIAL SCORE) 


* 12724; 


2 #* K 


* 12726; 


<2 ** K ) - 1 


*** 1272 0; 


2 <K - 1) 


* 12722; 


<2 ** < K - 1)) - 1 
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* 


12724 ) 


2 * * 


< K - 2) 


* 


12726 ; 


<2 ** 


G! ) - 1 


* 


12740 ; 


DCL + 


K - 1 


* 


12742 ; 


2 *.* 


< DCL - 1) 


* 


12744 ; 


S C 0 R E < I ) - S S E Q < I ) < SUBTRACT ADDRESSES) 


* 


12746 ; 


SSEQCJ) - SSEQCI) < SUBTRACT ADDRESSES) 


* 


12750; 


N X M 






12752; 


NEXT 


SCORE ADDRESS 




12754; 


NEXT 


SSEQ ADDRESS 




1 2 7 5 6 ; 


STORAGE LOCATION FOR A COUNTER VALUE OF 




12760; 


G < 1 ) < 


GENERATOR SEQUENCE 




12762; 


g ■: 2 ) , 


GENERATOR SEQUENCE 




12764; 


G < 2 ) , 


GENERATOR SEQUENCE 




12766; 


G < 4 ) » 


GENERATOR SEQUENCE 




12770; 


G < 5 ) » 


GENERATOR SEQUENCE 




12772; 


G < 6 ) , 


GENERATOR SEQUENCE 




12774; 


G < 7 >, 


GENERATOR SEQUENCE 




12776; 


G < 8 ) , 


GENERATOR SEQUENCE 



THE WORD COUNT IN LOCATION 1271 0 IS PLACED THERE BV 
RN INSTRUCTION RT THE END OF THE MESSAGE INPUT SUB- 
PROGRAM. THE VALUE OF THE GENERATOR SEQUENCES IS 
E N 1 E R E D B V 7 H E P R 0 G R f I M M E R B E F 0 R E S T A R T I N G R R UN. 
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APPENDIX B 



PROGRAM LISTING OF NOISE SIMULATION 

The following flow chart depicts the flow of the 
instructions in the following machine language listing 
of the noise simulation subprogram. The numbers (base 8) 
on the upper left of the blocks in the flow chart corres- 
pond to those instruction addresses of the subprogram 
mentioned in that block. 
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(11200-11324) 




This segment uses the noise para- 
meters placed in locations (10230) 
and (10236) to change the value 
of the number of set bits in a 16 
bit random number, so that the 
density of errors (XT) in the noise 
is varied. The noise sequence is 
stored in a block beginning at an 
| address in (13606) . ‘ 

(11550-11646) 
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The following computer printout is the noise program 
just discussed. 



011200 701670 0 
011202 7 0 0 2 3 7 6 
011204 7016701 
011206 7 0 0 2 3 7 4 
011210 7 0 1 6 7 0 2 
011212 7 0 0 2 5 2 4 
011214 7022767 
011216 ✓000001 
0112 2 0 7 0 0 2 4 6 2 
01122 2 / 0 0 1 0 O 2 
011224 ✓012021 
011226 7 0 7 7 2 0 2 
011230 ✓000426 
011222 ✓ 0 1 2 7 0 3 
011224 7O00020 
011226 7012704 
011240 700O020 
011242 7016705 
011244 7002426 
011246 7006210 
011250 7102405 
011252 7006211 
©11254 7005202 
011256 7001410 
011260 7077504 
0112 6 2 7 0 0 0 41 6 
011264 7006311 
011266 7005211 
011270 7005202 
011272 7001406 
011274 7077505 
01127 6 7 n 0 0 41 0 
0112 0 0 7 0 0 5 7 21 
011202 7012702 
0112 0 4 7 0 0 0 0 2 0 
011206 7000764 
011310 7005721 
011212 7012703 
011214 7O00020 
011216 7000766 
011220 7077420 
011222 7005720 
011224 7077224 
011326 7 0 0 0 0 0 0 



011330 


701670 0 


011232 


7002254 


011324 


7012701 


011226 


7004704 


011240 


7 012 7 0 2 


011342 


7 0 0 0 4 01 


011244 


7012727 


011246 


7044444 


011250 


7177204 


011352 


7 010 23 7 


011354 


✓177306 


011256 


7 0 0 5 2 2 7 


011360 


7177204 


011362 


7013720 


011364 


7177204 


011366 


7077107 


011370 


7 0 1 0 0 6 7 


011272 


✓002300 


011274 


7000240 


011376 


70 0 5 0 6 7 


011400 


7 0 0 2 2 6 4 


011402 


7 0 0 5 0 6 7 


011404 


7002262 


011406 


7016700 


011410 


7002176 


011412 


7016701 


011414 


7 0 0 217 0 


011416 


7 012 7 0 2 


011420 


7 0 0 0 0 2 0 


011422 


7012704 


011424 


7 0 0 0 0 2 0 


011426 


7 0 1 2 0 0 6 


011420 


7 0 2 0 0 6 7 


0114 2 2 


7 0 0 2 2 4 0 


011424 


7 0 0 1 0 0 2 


011426 


7 016 7 0 0 


011440 


7 0 0 214 6 


011442 


7012006 


011444 


70O5005 


011446 


7 0 0 6 2 0 6 


011450 


7005505 


011452 


7077402 


011454 


7005767 


011456 


7002224 
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011460 


7 0 01 0 0 5 


011574 


7006205 


011462 


701 6 704 


011576 


7006204 


G11464 


7002224 


011600 


7005505 


011466 


70062O5 


011602 


7 0 0 5 2 0 2 


011479 


7 0 7 7 4 0 2 


011604 


7001412 


011472 


7 0 0 0 4 0 7 


0116O6 


7077206 


011474 


7016727 


011610 


7 0 6 0 5 6 7 


011476 


7 0 0 2 2 1 2 


Oil 6 1 2 


7 0 0 2 O 5 2 


011500 


7177204 


011614 


7016702 


011502 


7010527 


011616 


7 0 0 2116 


011504 


7177206 


011620 


7006202 


0115O6 


7012705 


011622 


7020502 


01151O 


7177204 


011624 


7 0 0 2 7 6 0 


011512 


7005202 


011626 


7005267 


011514 


700240 5 


011620 


7002026 


011516 


7006211 


011622 


7 0 0 0 7 5 5 


011520 


7005205 


011624 


7 012 0 0 4 


611522 


7002272 


011626 


7012702 


011524 


7005211 


011640 


7000020 


011526 


7000725 


011642 


7077117 


011520 


7005721 


011644 


7600240 


811522 


7020167 


011646 


7 0 0 0 2 4 0 


011524 


7002052 


011650 


7016700 


011526 


7001001 


011652 


7001722 


O11540 


7 0 0 0 4 0 2 


011654 


7016701 


011542 


7012702 


011656 


7001724 


011544 


7000017 


011660 


7 0 1 6 7 0 2 


0115 4 6 


7 0 0 0 7 6 2 


0116 6 2 


7 0 017 2 4 


011550 


7 016 7 0 0 


011664 


7016705 


011552 


7 0 0 2 0 2 2 


011 6 6 6 


7002020 


011554 


70167O 1 


01167O 


7 0 1 2 0 0 2 


011556 


7 0 0 214 0 


011672 


7 012104 


01156O 


7 O12 7 0 2 


011674 


70 7 4 2 0 4 


011562 


7 0 0 0 0 2 0 


011676 


7O10422 


0115 6 4 


7 012O04 


011700 


7077505 


0115 6 6 


7 016 7 0 2 


811702 


7 0 0 0 0 0 0 


011579 


7002112 


* 




011572 


7 0 0 5 0 0 5 
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APPENDIX C 



Program Listing of Rate 1/n 
Convolutional Encoder 

The following flow chart depicts the flow of the 
rate 1/n convolutional encoder implementation. The cor- 
responding machine language program follows the flow 
chart and the instruction addressed are indicated on top 
of the respective flow chart block to which they corres- 
pond in the subprogram. 



( 10710 - 10726 ) 
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The following computer printout is the encoder implementa 
tion just discussed. 



010710 7016700 
010 712 7 0 0 2 6 6 4 
010714 7016701 
010716 7002662 
010720 7012704 
010722 7000020 
010724 7016706 
010 7 2 6 70 0 2 7 0 0 
010726 7016716 
010 7 2 2 70 0 2 754 
010724 7012027 
010726 7177204 
O10740 7005027 
010742 7177202 
010744 7000406 
010 7 4 6 7 0 1 2 7 4 6 
010750 7177202 
010752 7012027 
010754 7177204 
010756 7012627 
010760 7177202 
010762 7012746 
010764 70O0020 
010766 7012727 
010 7 7 0 7 0 0 0 0 0 1 
010772 7177214 
010774 7015746 
010776 7002700 
0110 0 0 7 012 7 0 5 
011002 7177202 
011004 7015702 
0110 0 6 7 0 0 2 6 1 6 
011010 7005105 



011012 


7 012 2 0 2 


011014 


7040502 


011016 


7 0 0 5 0 4 6 


0110 2 0 


7016746 


0110 2 2 


7 0 0 2 6 5 6 


011024 


7 0 0 5 7 2 6 


011026 


7 0 0 6 2 0 2 


0110 2 0 


700551 6 


011022 


7005246 


011024 


7 0 012 7 2 


011026 


7 0 0 6 211 


011040 


7005726 


011042 


7006226 


011044 


7005511 


011046 


7 0 0 5 20 4 


011050 


7 0 0 1 0 0 2 


011052 


7012704 


011054 


7 0 0 0 0 2 0 


011056 


7005721 


011060 


7005216 


011062 


7001252 


011064 


7005726 


011066 


7 0 0 5 216 


011070 


7 0 012 2 6 


011072 


7005726 


011074 


7005216 


011076 


7 0 0122 2 


011100 


7005704 


011102 


7 0 6 .1 4 O 2 


011104 


7 0 06 211 


01110 6 


7 0 0 5 2 0 4 


011110 


7001275 


011112 


7 0 0 0 0 0 0 
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APPENDIX D 



Program Listing of Viterbi Decoder 

The following flow chart depicts the flow of the 
Viterbi decoding algorithm implementation. The correspond- 
ing machine language program follows the flow chart and the 
instruction addresses are indicated on top of the respec- 
tive flow chart block to which they correspond in the 
subprogram. 
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( 11710 - 12042 ) 




( 12370 - 12546 ) ^ 




(*) 
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The following computer printout is the decoder implementa- 
tion just discussed. 



U11710 - 


7016700 


012050 


7016701 


011712 


7002014 


012052 


7001650 


011714 


7016701 


012054 


7 016 7O2 


011716 


7001706 


012056 


7001550 


011720 


7016702 


012060 


7 016 7 0 2 


011722 


7 0 01 6 6 6 


012062 


7 0 0 1 6 1 4 


011724 


7 016 70 4 


012064 


7 0 1 6 7 0 4 


011726 


7 0 017 0 0 


0 1 2 0 6 6 


70O1614 


01172O 


7 0 0 5 0 0 2 


012070 


700501 2 


011722 


7005012 


012072 


7 0 0 6 2 1 0 


011724 


7016714 


012074 


710 2 4 0 2 


011726 


7001740 


012076 


7 0 0 6 212 


011740 


7005102 


012100 


7077402 


011742 


7012105 


012102 


7 0 0 04 02 


011744 


7 0 4 0 2 0 5 


0121O4 


7006212 


011746 


7016744 


012106 


7005212 


011750 


7001720 


012110 


7077402 


011752 


7 0 0 5 0 0 6 


012112 


7 0 0 5 7 4 2 


011754 


7006205 


012114 


7077215 


011756 


7005506 


012116 


7005722 


011760 


7005 2 1 4 


012120 


7 016 7 0 2 


011762 


7001274 


012122 


7 0 0 155 4 


011764 


7 0 O 6 2 1 2 


012124 


7005010 


011766 


7006206 


012126 


7016704 


011770 


7005512 


012120 


7001552 


011772 


7005724 


0121 2 2 


7006210 


011774 


7005214 


012124 


7077402 


011776 


7 0 0 12 61 


012126 


7 0 6 221 0 


012000 


7 016 7 01 


012140 


7 0 7 7 2 0 6 


012002 


700162 2 


012142 


7 0 0 5 7 2 0 


012004 


7 O 0 5 72 2 


012144 


7005742 


012006 


7 0 0 5 1 0 2 


012146 


7077126 


012010 


7 0 2 6 7 0 2 


012150 


701670O 


0 1201 2 


7 0 0 1 7 1 2 


012152 


7 0 014 2 6 


012014 


7001407 


012154 


7016701 


0 1201 6 


7 0 2 O 2 0 0 


012156 


7 0 014 2 4 


0 1 2 0 2 0 


7 0 0 2 0 0 2 


012160 


7016702 


0 1 2 0 2 2 


7 0 6 0 0 0 2 


012162 


7 0 014 22 


012024 


7 0 0 0 7 4 2 


012164 


70167O2 


0 1 2 O 2 6 


7 1 6 0 0 0 2 


0121 6 6 


7 0 014 2 0 


0 1 2 0 2 0 


7 0 0 5 2 0 2 


01217O 


7 U 0 5 0 21 


0 1 2 0 2 2 


700O727 


012172 


7005O22 


012024 


7022767 


012174 


7016704 


012026 


7000001 


012176 


7001522 


0 12 0 4 0 


7 O01642 


012200 


70O5021 


012042 


7O01442 


0122O2 


7016722 


012044 


701670O 


0122O4 


7001514 


0 12 0 4 6 


7001542 


012206 


7O77404 
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012210 


/ 0 1 6 7 0 1 


012250 


7065267 


012212 


700140 0 


0122 5 2 


7601272 


012214 


7016702 


012254 


7001612 


012216 


7001276 


012256 


7005724 


O122 2 0 


7 0 1 6 7 0 4 


612260 


7005267 


01222 2 


7 0 012 5 6 


012262 


7 0 012 5 6 


612224 


7016767 


0 1 2 2 6 4 


7 0 0 1 0 0 2 


01222 6 


700146 0 


012266 


7 0 0 0 0 0 0 


012220 


7001410 


012270 


7 0 0 O 2 4 0 


6 1 2 2 2 2 


7 0 1 2 7 6 7 


012272 


7 0 0 0 2 4 0 


012224 


7 0 0 0 0 2 0 


012274 


7 0 1 2 7 6 7 


012 22 6 


7001404 


012276 


7 0 0 0 0 2 0 


012240 


7012227 


012400 


7081244 


012242 


7 1 7 7 2 0 4 


012402 


7010102 


012244 


701276 7 


012404 


7066702 


012246 


7 0 0 0 0 2 0 


012406 


7 0 012 2 4 


012250 


7 0 0 1 2 7 4 


012410 


7 0 0 5 0 27 


012252 


7005267 


012412 


7177202 


012254 


700146 2 


012414 


7 016 7 0 6 


012256 


7 0 0 20 51 


012416 


7001276 


012260 


7016767 


012420 


7012727 


012262 


7001444 


012422 


7 0 0 0 0 0 1 


012264 


7001262 


012424 


7 1 7 7 2 1 4 


012266 


7 0 1 0 2 6 7 


012426 


7005267 


012270 


7001260 


012420 


700121 2 


01227 2 


7012267 


0124 2 2 


7801402 


012274 


7001256 


012424 


7077607 


012276 


7005267 


012426 


7000414 


012200 


7001246 


012440 


7 012 7 6 7 


012202 


7001405 


012442 


7 0 0 0 0 2 0 


012204 


7021267 


012444 


7601176 


0122O6 


7001244 


0124 4 6 


7012767 


012210 


7 0 0 2 7 6 6 


012450 


7 1 7 7 2 0 2 


012212 


7 0 0 5 7 2 2 


012452 


760115 6 


012214 


7 0 00 77 0 


012454 


7012227 


012216 


7166767 


012456 


7177204 


0 1 2 2 2 0 


7 0 014 2 2 


012460 


7 016 7 27 


012222 


7 0 012 2 6 


012462 


7 0 6114 6 


012224 


7 016 70 6 


012464 


7177202 


012 2 2 6 


7 U 0 12 2 2 


012466 


7 0 0 0 7 6 2 


0122 2 0 


7 0 0 6 2 1 4 


012478 


7 0 1 6 7 6 7 


01222 2 


7026716 


012472 


7 O 012 4 0 


012224 


7 0 014 0 4 


012474 


7001160 


012226 


710 0 0 0 4 


012476 


7 G 1 0 2 6 7 


0 1224 0 


70 0 0 24 0 


012500 


7001250 


012242 


700O240 


012502 


7 0 1 0 1 6 7 


012244 


7000240 


012504 


7001246 


012246 


7005214 


0125O6 


7O20267 
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012516 


7001106 


012650 


7000240 


012512 


7100010 


012652 


7066767 


0125-14 


7066767 


012654 


7000756 


012516 


7 0 0 1 2 2 6 


012656 


7001000 


012520 


7001250 


012660 


7010267 


012522 


7066767 


012662 


7000746 


012524 


7001 2 2 0 


012664 


7066767 


012526 


7 0 01 2 2 4 


012666 


7001040 


012550 


7 0 0 0 4 0 7 


O12670 


7 0 00 7 4 0 


012552 


7000240 


012672 


7017767 


012554 


7166767 


012674 


7 0 U 0 7 5 4 


012556 


7 0 01206 


012676 


7000762 


012540 


700121 0 


012700 


7016706 


012542 


7 1 6 6 7 6 7 


012702 


7000774 


012544 


7001200 


012704 


7015767 


012546 


7001 2 0 4 


O12706 


7177502 


01255 0 


7 0127 6 7 


012710 


7000726 


012552 


7 0 0 0 0 0 2 


012712 


7012005 


012554 


7 0012 0 0 


012714 


7074567 


012556 


7011267 


012716 


7 0 0 6 7 2 0 


012560 


7001076 


012720 


7 0 0 5 0 6 7 


012562 


7 016 7 0 6 


012722 


700071 0 


012564 


7001112 


012724 


7005167 


012566 


7015767 


012726 


7000710 


012570 


7177502 


012750 


7016705 


012572 


7 00104 4 


012752 


7000750 


012574 


7012005 


012754 


7016767 


012576 


7074567 


0127 5 6 


7000776 


0126 00 


7001056 


012740 


7000674 


O12602 


70O5067 


012742 


7046767 


012604 


7 0 0 1026 


012744 


7000672 


012606 


7 0 0 516 7 


012746 


7 0 0 0 6 6 6 


012610 


7001 0 2 6 


012750 


7066767 


0 1261 2 


7 0 1 6 7 0 5 


012752 


7 0 0 0 6 6 2 


012614 


7001 0 6 6 


012754 


7000656 


012616 


7016767 


012756 


7 0 O 6 2 6 7 


012620 


7001114 


O12760 


7 0 0 0 6 5 6 


0 1 2 6 2 2 


7 0 0101 2 


012762 


7 0 7 7 5 0 5 


0 1 2 6 2 4 


7 0 4 6 7 6 7 


012764 


7 07 7617 


012626 


7001010 


0 1 2 7 6 6 


7066767 


012650 


7 0 6 10 0 4 


O12770 


7 0 O064 2 


012652 


7066767 


012772 


7000666 


012654 


7 0 O 10 0 0 


012774 


7 0 2 6 7 6 7 


012656 


7 0 0O774 


012776 


7 O 0 0 6 6 0 


012640 


7 0 0 6 2 6 7 


015 0 0 0 


7 O 0 0 6 6 0 


012642 


7 0 0 0 7 7 4 


0 1 5 0 0 2 


7005421 


012644 


7077505 


015004 


7016705 


012646 


7077617 


015006 


7000742 
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0 1 :< 0 i e 


7016715 


012100 


7010667 


012012 


7000646 


0121O2 


7 0 0 0 6 5 0 


012014 


7010167 


012104 


7O05 2 6 7 


012016 


7 0 0 0 6 4 4 


012106 


7 O 0 0 6 4 6 


012020 


7066767 


012110 


7 0 0 1 2 2 2 


012022 


7 0 0 0 7 0 4 


012112 


7022122 


012024 


7 0 0 0 6 2 6 


012114 


7000240 


012026 


7 0 1 6 7 0 5 


012116 


7005267 


012020 


7 0 0 0 6 2 2 


012120 


7000524 


012O22 


7016706 


012122 


7O01212 


0 1 2 0 2 4 


7 0 0O716 


012124 


70167O0 


012O26 


7 011516 


012126 


700D462 


01204O 


7 0 0 6 2 1 6 


012120 


7 0 2 0 2 6 7 


012042 


7005216 


0121 2 2 


7000464 


012044 


700041 0 


012124 


7 100406 


O12046 


7016705 


012126 


7016701 


012050 


7 0 0 0 7 0 0 


012140 


7000452 


012052 


7 016 7 0 6 


012142 


7 O16 7 0 2 


012054 


7 0 0 0 6 7 6 


012144 


70O0450 


012056 


7016715 


012146 


7 1 6 2 7 0 7 


012060 


700O576 


012150 


7000700 


012062 


7011116 


012152 


7016701 


012064 


7006216 


012154 


7000442 


012O66 


7 0 1 6 7 0 5 


012156 


7016702 


012070 


7 0 0 0 6 6 0 


012160 


7 0 6 0 4 4 0 


012O72 


7 0 2 2 5 2 6 


012162 


7 162707 


012074 


7 0 1 0 5 6 7 


012164 


7 0 0 0 714 


0 12 0 7 6 


7 0 0 0 6 5 2 


012166 


7 0 0 0 0 0 0 
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APPENDIX E 



Program Listing of Supplementary Subprograms 

1. Initialization subprogram (10174-10706) 

This subprogram enters the necessary addresses and 
constants into those locations noted in Appendix A by * . 

The computer listing below is the initialization subprogram. 



O10174 


7010700 


010176 


7 0 1 2767 


010200 


7000002 


010202 


7003474 


610204 


7 0 1 2 7 6 7 


0102 0 6 


7 0 0 0 0 0 3 


010210 


7003470 


010212 


7612767 


010214 


7000001 


010216 


7003464 


010220 


7 0 1 2 7 6 7 


610222 


7 0 0 0 0 2 0 


010224 


7 0 0 3 4 6 0 


010226 


7 012767 


010230 


7 0 0 0 0 0 1 


010232 


7 003456 


010234 


7012767 


010236 


7 0 0 0 0 0 0 


010240 


7 0 0 3 4 5 2 


010242 


701O701 


010244 


7062701 


0 10 24 6 


7 003514 


010250 


7 0 1 0 1 6 7 


0 10 2 5 2 


7 0 0 2 ± 5 2 


010254 


7 0 6 2 7 0 1 


010256 


7 0 0 0 4 0 6 


010260 


7010167 


010262 


7 0 U i ^ 4 4 


010264 


7 0 6 2 7 0 0 


010266 


7 0 0 3 4 0 2 



010270 


7010 


001 


010272 


70 6 2 


70 0 


010274 


7000 


6 0 0 


010276 


7010021 


010300 


7 0 6 6 


700 


0 1 © 3 0 2 


70 0 3 


4 04 


010304 


7010 


021 


01O306 


7 016 


73 7 


010310 


7003 


3 6 6 


010312 


717 7 


304 


010314 


7016 


737 


010316 


7 0 0 3 


3 70 


010320 


7177 


306 


010322 


7013 


76 7 


010324 


7177 


30 4 


010326 


70 03 


4 20 


010330 


706 3 


700 


010332 


7177 


30 4 


010334 


7010 


021 


010336 


7016 


73 7 


0 10 3 4 0 


7 0 0 3 


342 


010342 


717 7 


3 0 6 


010344 


7 0 1 3 


767 


0 1034 6 


717 7 


3 0 4 


010350 


7003 


34 6 


010352 


70 6 3 


700 


010354 


7177 


30 4 


010356 


7010 


0 21 


6 1 0 3 6 0 


7 0 6 3 


7 00 


010362 


7177 


30 4 
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010264 


7010021 


010526 


7016702 


010266 


7062700 


010540 


7002142 


010270 


/I 7730 4 


010542 


7 0 0 5 0 0 2 


010272 


/ 0 1 0 0 2 1 


610544 


7 0 0 5 2 0 2 


010274 


7 0 0 5 0 0 2 


010546 


7 0 0 6 2 0 2 


010276 


7 0 0 5 2 0 2 


010550 


7 0 7 7 2 0 2 


010400 


7016702 


010552 


7005202 


010402 


7 0 0 2 2 7 6 


010554 


7010267 


010404 


7 0 O 6 2 0 2 


010556 


7 0 0 215 6 


010406 


7 0 7 7 2 0 2 


0105 6 0 


7 0 0 5 2 0 2 


010410 


7010 2 6 7 


010562 


7006202 


010412 


7002210 


010564 


7010267 


010414 


7 O 0 6 2 0 2 


010566 


7002106 


010416 


7 0 6 0 2 0 0 


010570 


7 0 1 6 7 0 2 


010 4 2 0 


701002 1 


G10572 


7 0 0 2 1 0 6 


010422 


7006202 


010574 


7 0 6 5 2 0 2 


010424 


7 0 0 5 2 0 2 


010576 


7066702 


010426 


7010267 


01060O 


7002104 


010420 


7002274 


010602 


7010267 


010422 


7005202 


010604 


7 0 0 2 1 2 2 


010424 


7 012 7 0 2 


010606 


7016702 


010426 


7 0 0 0 0 0 2 


010610 


7 0 0 20 74 


010440 


7060200 


010612 


7005202 


010442 


7010021 


010614 


7 0 0 5 0 0 2 


010444 


7077202 


010616 


7 0 6i 5 2 O 2 


010446 


7 0 0 6 2 0 2 


010620 


7 0 0 6 2 0 2 


010450 


7 010 2 6 7 


010622 


7077202 


010452 


7002254 


010624 


7010267 


010454 


7005202 


0 1 0 6 2 6 


7 0 0 211 2 


010456 


7010267 


010620 


7016702 


010460 


7002250 


0 1 0 6 2 2 


7 0 0 2 7 6 4 


010462 


70 0 5 2 0 2 


010624 


7166702 


010464 


7 0 0 6 2 0 2 


0 1 0 6 2 6 


7 0 0 2 7 5 6 


010 4 6 6 


7 010 2 6 7 


0 1 0 6 4 0 


7010267 


010470 


7 0 0 2 2 4 2 


010642 


7 0 0 21 0 0 


010472 


7 0 1 2 7 2 7 


010644 


7 016 7 0 2 


01O474 


7 0 0 0 O 0 4 


0 10 6 4 6 


7 0 6 27 5 2 


010476 


7177204 


0 1 0 6 5 0 


7 1 6 6 7 O 2 


0 1 0 5 0 0 


7016727 


010652 


7 9 0 2 7 4 2 


010502 


7 0 a 2 17 6 


010654 


7 010 2 6 7 


010504 


7 1 7 7 2 0 6 


0 1 0 6 5 6 


7 0 0 2 0 6 6 


010 5 0 6 


7012767 


0 1 0 6 6 0 


7 010 7 0 0 


0 1051 0 


7177204 


0 1 0 6 6 2 


7 0 6 2 7 0 0 


0 1051 2 


70022O6 


0 1 0 6 6 4 


7 0O1 2 6 4 


010514 


7016727 


010 6 6 6 


7 010 O 6 7 


010516 


7002164 


010670 


7002204 


0 1052 U 


7177204 


010672 


7 0 0 6 26 7 


010522 


7016727 


010674 


7002012 


010524 


7002152 


O10676 


7006267 


010526 


7 1 7 7 2 O 6 


010 7 0 0 


7O02016 


010520 


7 0 1 2 7 6 7 


010702 


7 0O6267 


010 5 2 2 


7177204 


0107O4 


7002042 


010 5 2 4 


7002160 


010706 


7 0 O 0 0 0 0 
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2. Message Input Subprogram (10000-10144) 

This subprogram was used to indicate that an alphanu- 
meric symbol typed at the Keyboard was entered into core 
(echo) . The symbols are stored in core in an ASC II 
(7-bit) code representation. The Keyboard symbol, @ , 
is used to terminate message entry and the number of 
8-bit computer bytes used is stored in location (13710 
Appendix A) . The following computer printout is a listing 
of the ASC II message input subprogram. 



e i e 0 o £i 


2010700 ' 


010064 


7 0 £1 0 


012 


0 1 £i £i 0 2 


7062700 


010066 


7105 


7 27 


£i 1 £i £i 0 4 


7004376 


010070 


7177 


56 4 


61 £100 6 


7005002 


010072 


710 0 


27 5 


010010 


7105737 


010074 


7112 


727 


0 1061 2 


7177560 


010076 


7 0 0 0 


20 0 


010014 


7100375 


010100 


7 1 7 7 


56 6 


010016 


7113710 


010102 


7077 


10 7 


010020 


7 17756 2 


010104 


7105 


7 2 7 


010022 


7122710 


010106 


717 7 


56 4 


010024 


7 0 0 0 3 0 0 


010110 


7 1 0 0 


27 5 


010026 


70 014 3 5 


010112 


7112 


72 7 


0 1 0 0 :< 0 


7105727 


010114 


7 0 0 0 


212 


0 1 0 0 3 2 


7177564 


010116 


717 7 


5 6 6 


0 1 0 0 :< 4 


7160375 


0101 2 0 


7 0 0 0 


732 


0 1 0 0 2 6 


7112 0 2 7 


010122 


7 0 1 2 


7 02 


0 10 0 4 0 


7177566 


010124 


7 0 0 0 


0 1 0 


010042 


7 0 0 5 2 £i 2 


010126 


710 5 


020 


010044 


7 123727 


0101 3 8 


7 0 6 5 


2 £i 2 


0 10 0 4 6 


7177562 


010132 


7077 


2 0 2 


0 1 0 0 5 0 


7 0 0 0 215 


010134 


7 0 0 6 


2 0 2 


0 1 0 0 5 2 


70O1256 


010126 


7 0 0 6 


20 2 


010054 


7112740 


010140 


7 0 1 0 


cl t* l 


010056 


7 0 0 0 2 4 £i 


010142 


7 0 0 2 


544 


010060 


7 1 £i 5 7 2 0 


010144 


70 0 0 


000 


010062 


7012701 


* 
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3. Message Output Subprogram (13200-13364) 

This subprogram is basically the same as that in 
section 2 of this appendix. However, besides typing what 
is in core starting at a location specified at (13602) , 
the number of errors (bit differences) are determined 
between the input message and the decoded out message. 

This message error value is stored at (13672) (Appendix A). 
The following computer printout is a listing of the ASC II 
(7-bit) code output subprogram. 



012200 


2016700 


012274 


7100275 


013202 


7 0 0 027 6 


013276 


7112727 


013204 


/ 012701 


013200 


700001 2 


012206 


/ 0 0 0 0 4 0 


012302 


7177566 


012210 


2111027 


013204 


7105737 


013212 


2177566 


013206 


7177564 


012214 


2122720 


012310 


7100275 


012216 


2000000 


013212 


7000734 


G12220 


2001425 


013314 


700O240 


013222 


2 105737 


013216 


7000240 


012224 


2177564 


013320 


7 016 7 0 0 


012226 


2100375 


012322 


700O254 


012230 


2 0 0 5 301 


013324 


7016701 


013232 


210 0 2 6 6 


012326 


7 0 0 6 2 5 2 


013234 


2122710 


013330 


7 0 0 5 0 0 2 


013236 


2000240 


013332 


7 011 0 0 2 


O13240 


2 0 0 1 3 6 2 


012334 


7 0 7 4 211 


013242 


2 1 1 2 7 2 7 


013326 


7012704 


013244 


2 0 0 0 015 


012240 


7000020 


013246 


217 7 5 6 6 


013242 


7 0 0 6 211 


012250 


2 6127 0 2 


012244 


7 0 0 5 5 0 2 


013 2 5 2 


2 0 0 0 012 


013246 


7 0 7 7 4 0 2 


013254 


2105727 


0132 5 0 


7 0 2 2 0 2 1 


012256 


2 1 7 7 5 6 4 


013252 


7020067 


013260 


2 100275 


012254 


7000224 


013262 


2 1 0 5 0 2 7 


012256 


7001365 


013264 


2177566 


01326O 


7010267 


013266 


2077206 


012262 


7000206 


013270 


2105737 


012264 


7 0 0 0 0 0 0 


013272 


2 1 7 7 5 6 4 


+ 
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4. Analysis Subprograms (13250-13544) 

From address (13250) to address (13376) a short sub- 
program is listed that is used to step through well defined 
encoders of rate 1/2, 1/3, and 1/4. In the last addresses 
(13400-13544) a program is listed which determines the 
number of error bits in 16-quantized bits and stores the 
results so that a distribution of errors may be plotted 
as in figures 17-23. (See next page.) 
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012400 


7016700 






012402 


7 0 0 ©202 






012404 


7016701 






012406 


7O00226 


012250 


7005267 


012410 


7006201 


012252 


700051 2 


012412 


70O0240 


012254 


7 0 2 6 7 2 7 


012414 


7 0 1 2 7 6 7 


012256 


7000506 


012416 


7 0 0 G 0 2 0 


012260 


7 0 0 0 0 1 0 


012420 


7000214 


012262 


7001 0 2 0 


G12422 


7012704 


012264 


7012767 


012424 


7014100 


012266 


7 0 O 0 0 0 4 


012426 


7012705 


012270 


700O474 


012420 


7 0 0 0 0 2 0 


012272 


700 5 2 6 7 


012422 


7 0 0 0 2 4 0 


012274 


7 0 0 0 4 6 6 


012424 


7000240 


012276 


7 0 2 6 7 2 7 


012426 


7 0 0 0 2 4 0 


O12200 


7 0 d 0 4 6 2 


012440 


7 0 0 0 2 4 0 


012202 


7 0 0 0 0 1 0 


012442 


7 0 0 5 0 0 6 


012204 


7001017 


012444 


7G16702 


012206 


7 0 1 2 7 6 7 


©12446 


7000224 


012210 


7 0 0 0 0 0 4 


012450 


7 0 0 5 0 0 2 


012212 


7000450 


012452 


70O6202 


012214 


7005267 


612454 


7 0 0 6 2 1 0 


012216 


7000442 


012456 


700550 2 


012220 


7026727 


O12460 


7005267 


01222 2 


700O426 


012462 


7000152 


012224 


700O020 


012464 


7001407 


012226 


70 01 0 0 6 


012466 


7O77207 


0122 2 0 


7012767 

| 


012470 


70202O1 


012222 


7 0 0 0 0 1 0 


012472 


7002401 


012224 


7000424 


D12474 


7005206 


0 1222 6 


7062767 


012476 


7005205 


012240 


7 0 0 0 0 0 2 


012500 


7001414 


012242 


7000414 


012502 


7000760 


012244 


7016777 


012504 


70O5720 


012246 


7 0 0 0 2 2 2 


012506 


7026700 


01225 0 


7 O 6 0 41 6 


0 1251 0 


7000076 


012252 


7 0 6 2 7 6 7 


0 1251 2 


70O2002 


012254 


7 0 0 0 0 O 2 


012514 


7 0 G 0 0 G 0 


012256 


7 0 0 0 410 


D12516 


7 6 O G 2 4 0 


0122 6 6 


7 O 2 6 7 2 7 


012520 


7000240 


012262 


7 8 0 0 2 7 4 


012522 


7012767 


012264 


7 O 0 0 0 2 1 


012524 


7 0 0 0 0 2 0 


0 1 2 2 6 6 


7 0 014O2 


0 1 2 5 2 6 


7000106 


012270 


7 0 0 012 7 


0 12 52 0 


7000756 


CH2272 


7010000 


012522 


70062O6 


012274 


70O0127 


012524 


7060604 


01227b 


70D1172 


612526 


7005214 






012540 


7O00240 






012542 


700O727 






012544 


7O00000 
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APPENDIX F 



SAMPLE RUN 



THE COMPUTER LISTING BELOW THIS PARAGRAPH IS A 
DEMONSTRATION OF THE PROCEDURE REQUIRED TO USE THE 
ENTIRE PROGRAM PRGPERLV. THE PROGRAM IS STORED ON 
A DISK UNDER THE NAME VSAC. SAV. THE CONVOLUTIONAL 
CODE USED IS A RATE 1X2- (l< = 3> CODE. OTHER INPUTS 
ARE G=l> DCL=16 (20 IN BASE 8>» AND THE GENERATOR 
SEQUENCES OF 3 AND 4 (BASE &>. 



~C: 

. GET VSAC. SAV 

. START 1172 
ODT V61-01 
+ 

*10200/000002 

*10206/000003 

*10214/000001 

* 10222/000020 

*10230/000001 

+10236/000001 

+ 

*12336/100004 

* 

+ 1376 0 / 0 0 0 0 0 3 
+ 13762 / O 0 0 0 O 4 
* 

* ~c 

STftRT 10000 
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The paragraph below is the input message for the sample 

run . 



THIS IS Ft SAMPLE RUM OF A THESIS COMPUTER 
PROGRAM USED TO EVALUATE THE PERFORMANCE OF RATE 
IV N CONVOLUTIONAL CODES, OVER A SIMULATED NOISV 
CHANNEL, WITH A VITERBI DECODER. THE NOISE PARAMETERS 
FOR THIS RUN (CORE LOCATIONS 10220 AND 10226) ARE 
THE SAME AS WERE USED TO OBTAIN THE DISTRIBUTION 
IN FIGURE 20 OF THE THESIS TEXT. THIS MESSAGE WILL 
NOW BE TERMINATED BV THE ' AT' CHARACTER ON THE PANEL. 



When the message above is encoded (as given above) and the 

noise is added, the resulting decoded message is: 

THIS IS! A"S0MPLE RUN OF E THESIS 
COHPUTR PROGBAM USED SO EVALUATE 
THE ERFNRUANCEDOF RATE 1/N CONVOLUTIONAL ! DCDES, 

0 V M R APIMULBTEE NOISV PANNEL ITH<0 
VITERBI DECOLOR. UHE NOISE PARAMETERS 
FOR TPIS UN (COU LOCATIONS 10220 
AND 102261 ARU TE SAUA AS TEVE USED "TO 
OBTAIN THE DISTRIBUTION I HOF JURE 
21 OF THE THEIP TEXD. THIS MESSAGD 
WILL NIW BE TERMINATED B V " T H M 'AX" 

CHBRCTER ON THM 0 B NED. 

When no coding is applied to the input message and the 
same noise is added the received message is: 

1L A P @ f< I N B 5 R D 8 B 0 J 0 M E 0 V U N " M D 0 0 < T I C K [ O C 0 M X U 
D R R 0 A M < U R M D t D _ < M T C M 3 A U E D D A 0 9 D R F 0 E A F K D " F 
S K U I " A ! R M V A V R M " E 6 C 0 E R > 0 D S D D ! N N A S D " T C R 0 I 
UUGRSONKV . <5— * 

PJIS ZUO ! HCOREfLBI T I 0J[ 1 0 : ? " E D " A 4 2 ; V A £ 

R & E B I 0 U f- 1 B f G R M B E C U D D T 0 0 M D V E I N ! D U D S T P H C W U I 
V0(?E@?B80JRU! BQTtifNOGvE T G R 
BJRI K P ft R " _NDJH ! XNEL. D 
* Improper carriage return received. 
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